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ABSTRACT 

We study the ability of three-frequency, arcminute-resolution microwave measurements to detect 
galaxy clusters via their Sunyaev-Zel'dovich (SZ) distortion of the microwave background. For this 
purpose, we have constructed large-area simulations of the microwave sky, and we have made them 
publicly available to further investigations into optimal data reduction techniques for upcoming SZ 
cluster surveys. In these sky simulations, galaxy clusters are modeled using N-body simulated dark 
matter halos plus a gas prescription for the intracluster medium that allows the small scale cluster 
physics such as star formation and feedback to be realistically incorporated. We also model the 
primary microwave background, radio and infrared point sources, galactic dust emission, and the SZ 
flux including kinetic and relativistic contributions. We make use of these simulations to study the 
scaling relation between integrated SZ flux and cluster mass, and find our clusters follow a power-law 
with an index that is steeper than that for self-similar cluster models. Some evolution of the power- 
law index and normalization with redshift is also observed. These simulations are also used to study 
cluster detection for the Atacama Cosmology Telescope (ACT). Using a multi- frequency Wiener filter 
to separate clusters from other microwave components, we find that ACT alone can recover a cluster 
sample that is « 90% complete and « 85% pure above 3 x 10 14 M Q . 

Subject headings: cosmology: theory — galaxies:clusters:general — intergalactic medium — mi- 
crowave : galaxies: clusters 



1. INTRODUCTION 

Cluster catalogs provide valuable information on the 
evolution and distribution of matter over cosmic time. It 
has long been realized that the abundance of galaxy clus- 
ters as a function of redshif t is a sensitive probe of the 
underlying cosm ology (e.g. lOukbir fc Blanchard [19921: 
Eke et all Il996l: IViana fc Liddld Il996t BarbosalEaD 
1996t iBahcall fe Fan! 119981 : iRosati et all l2002ft~ It 
has also been appreciated that a large sample of 
galaxy clusters can provide important constraints 
on the dark energy density and equation of state 
which arc complementary to those obtained from 
microwave ba ckground and type l a supernovae measure- 
ments (e.g. IWang fc Steinhardtl 119981: lHaiman et al.1 



200H: IWeher et all l2002t iWeller fc Battvd 120031: 



Maiumdar fc Mohrll2003t IWang et al.ll2004t iLima fcHul 
2005ft . This has prompted efforts to obtain cluster 
catalogs from wide, deep surveys. 

Searches for clusters based on their Sunyae v-Zel'dovich 
(SZ) effect (jSunvaev fc Zeldoviclil970l 1 1972ft have a par- 
ticular advantage over X-ray and optical searches since 
the SZ signal does not dim at high redshifts (z>l), where 



the cluster abundance is strongly de pendent on cosmol - 
ogy (For reviews on t he SZ effect see Birkins hawl ()1999f ) 
and ICarlstrom et al.l ([2002])). In light of this, there 
are a host of ground-based bolometer array instruments 
(ACBAR, ACT, APEX, SPT) and interferometers (AMI, 
AMiBA, SZA) either online or scheduled to come online 
in the next few years that potentially will detect hun- 
dreds to thousands of clusters at microwave frequencies 
via their SZ effect (iRunvan et alJ 120031: iKosowskv l 12003 : 
Giisten et aill2006l : feuhl et al.ll2004j : iKneissl et alJboOl : 
Li et al.ll2006ft . The Planck satellite, to be launched in 



2008, should also provide an all-sky map o f massive clus- 
ters in the microwave band (|Tauberll200l ). 

Using cluster catalogs to constrain cosmology requires 
a solid understanding of the cluster selection criteria. A 
given survey's cluster selection function is intimately tied 
to characteristics of the instrument and data reduction 
techniques. Considerable attention has been devoted to 
exploring various data reduction methods for different 
SZ surveys. This has been pursued in an effort to study 
both the optimal detection of clusters and the optimal 
recovery of the cluster SZ flux. The former consists of 
maximizing the completeness and minimizing the con- 
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tamination of a given cluster sample. The latter is nec- 
essary to properly identify the observable (SZ flux) with 
mass, or to employ self-calibration schemes to fit for the 
SZ flux - mass relation and cosmology simul taneously 
(jMaiumdar fc Mohdl200l I Younger et alJl2006f >. 

Various filtering approaches have been utilized to 
isolate the SZ signal and extract clusters in simula- 
tions. Such techniqu es have employed matched filters 
(iHerranz et al. 2002a, b ), Wiener filter s (Aghani m et al.l 
1997), wavelet filters ([Pierpaoli et al.l 120051). an d max- 
imum entropy methods (|Stolvarov et al.l l2002f ). Sev- 
eral papers have implemented versions of these filters 
and begun careful s tudy of selection functions for up- 
coming SZ surveys (ISchulz fc White! 120031: Iwliite] 1 20031: 



Melin et al.ll2005t iVale fc Whitell200a iMelin et al.ll200a 
Schaefer fc Bartelmannl l2006t Uuin et all I2007D . With 



this in mind, we have created a large-scale simulation of 
the microwave sky on which these data reduction tech- 
niques can be refined. 

To model the SZ flux accurately generally requires ex- 
pensive, high-resolution hydrodynamic simulations to re- 
alistically model the small scale cluster physics. How- 
ever, it is a challenge to create such simulations in a large 
enough volume so that one is not limited by cosmic vari- 
ance when performing statistical studies. To overcome 
this challenge, we create a cluster catalog using an N- 
bo dy simulation combin ed with a gas prescription given 
bv lOstriker et all (|2005fl . The N-body code only needs 
to be run once for a given cosmology, and the small scale 
cluster physics (including non-spherically symmetric gas 
distributions, star formation, and feedback) are added 
afterward via the gas prescription. As a result, the clus- 
ter physics can be varied easily without having to redo 
expensive runs. 

These cluster simulations were incorporated within 
microwave-sky simulations consisting of two strips of the 
sky, each 4 degrees wide in declination and 360 degrees 
around in right ascension. One strip is centered at a dec- 
lination of -5 degrees, and the other is centered at a decli- 
nation of -55 degrees, to match the portion of sky that the 
Atacama Cosmology Telescope (ACT) will observe. The 
galaxy clusters within this simulation have halo masses 
down to 5 x 10 13 M , and the N-body simulation employs 
the latest cosmological parameters derived from a combi- 
nation of WMAP3, SDSS, HST, and SN Astier observa- 
tions (fi m =0 26, n A =0.74 fi & =0.044, n=0.95, cr 8 =0.77, 
h=0.72) (see lSpergel et al.ll2006L and refer ences therein) . 
We include the full SZ effect, with relativistic corrections, 
as well as infrared and radio point sources ( uncorrelated 
with c lusters) , using number coun t s give n by lBorvs et al.l 
(|2003h (infrared) and lKnox etaD (|2004f ) (radio). Galac- 
tic dust and primary microwave background fluctuations 
are also m odeled, the later generated using the WMAP 
ILC map ([Bennett et al.l l2003al iHinshaw et ail l2006f ). 
These sky simulations were made at the observing fre- 
quencies of 145, 215, and 280 G Hz and are in a cylin- 
drica l equal area sky projection (jCalabretta fc Greisenl 
12001 ). The final sky maps are in a standard FITS for- 
mat and have a pixel size of « 0.2' x 0.2' and units of 
Jy/ster. A catalog of all the cluster halos in the simula- 
tion is also provided, specifying each halo's basic prop- 
erties. These microwave-sky simulations are available at 
|http:/ /www. a stro. princeto n.edu/~act/| 

We use these simulations to study the scaling relation 



between SZ flux and cluster mass, as well as prospects 
for cluster detection with ACT. Roughly 10 5 clusters are 
used to determine the simulated SZ flux - mass scaling re- 
lation. We compare cluster M200 to Y200, which we take 
to be the SZ Compton-y parameter integrated over a disk 
of radius R200 ■ These clusters are fit to a power-law re- 
lation between Y200 x E(z)~ 2 / 3 and M200, and estimates 
of the power-law index and normalization are given. 

We also employ a multi-frequency Wiener filter and 
simple peak-finding algorithm to forecast cluster detec- 
tion given ACT instrument specifications. Since there 
has been considerable interest of late in how point source 
contam ination affects cluster detection (e.g. I Melin et al.l 
(2006)), this is investigated under three different con- 
tamination assumptions: no point sources, only infrared 
point sources, and both infrared and radio point sources. 
Completeness of our recovered cluster sample is given as 
both a function of M200 and Y200 (in units of arcmin 2 ), 
the latter being directly obtainable from microwave ob- 
servations. We also give the purity of our projected clus- 
ter sample, where purity is one minus the percentage of 
false-positive detections. The issue of optimal SZ flux 
recovery for individual clusters is left to later work. 

The outline of this paper is as follows. In §2, we give 
the simulation details. §3 consists of an investigation of 
the Y-M relation as suggested by our simulations, and §4 
gives cluster detection projections for ACT. In §5, we dis- 
cuss directions for future work, and, in §6, we summarize 
and conclude. 

2. SIMULATIONS 

2.1. Simulated Clusters 

2.1.1. The Dark Matter Run 

Initial conditions for the TV-body run were generated 
with the GRAFIC2 code (|Bertschingerll2001h . available 
at Ettp : / /web. mit.edu/edbert/ As written, this 
code uses a spherical Hanning filter on small scales. 
However, we have found that this filter significantly 
suppresses power on these scales, so we removed it from 
use. GRAFIC2 perturbs particles from a regular grid 
using the Zel'dovich approximation. The simulation 
started at redshift z=35.3, when the initial density 
fluctuation amplitude on the scale of this grid was 
10%. N = 1024 3 w 10 9 particles were contained 
in a periodic box of size L = 1000/i _1 Mpc, so the 
particle mass is thus 6.72 x lO lo /i _1 M . The cubic 
spline softening length was set to e = 16.276/i _1 kpc. 
The simulat i on wa s carried out with the T PM code 
(|Bode et al.l 120001: iBode fc Ostrikerl 120031 ). modi- 
fied sligh tly from the publicly available version (at 
http : //www . astro . princeton . edu/$\sim$bode/TPM/) . 
Particle positions and velocities were followed at dou- 
ble precision, though acceleration was kept at single 
precision. Also, no lower limit was set to the param- 
eter B used in domai n decomposition (see Eqn. 5 of 
IBode fc Ostrikerl [2003) , which at late times results in 
more particles being followed at full force resolution. 
The initial domain decomposition parameters used are 
A = 1.9 and B = 9.2. The PM mesh contained 2048 3 
cells, and the maximum sub-box was 256 cells. By the 
end of the run, 5 x 10 6 trees containing 54% of the 
particles were followe d at full res o lution . More details of 
the simulation are in IBode et all ((20061) . 
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Fig. 1. — The halo abundance above a minimum mass per unit 
redshift per 100 square degrees measured from the simulation is 
compa red with the semi-analytic fitting formula from Jenki ns" et alj 
(2001). From top to bottom, the three sets of points are for the 
minimum masses Af m i n = (1,2,4) X 10 14 Mq /h. 



At each PM time step, the matter distribution in a thin 
shell is saved. The radius of the shell corresponds to the 
light travel time from az = observer sitting at the ori- 
gin of the box, and its width corresponds to the time step 
interval. The portion of a spherical shell covering one 
octant of the sky was saved, so for comoving distances 
larger than 1000/i _1 Mpc there can be some duplication 
in structures as the periodic box is repeated. However, 
while there are repeated dark matter halos in the simu- 
lations, the repetition is usually at a different redshift, so 
the dynamical state is different. In cases where a cluster 
appears twice at the same redshift, it is viewed at two 
different angles, making each projected image unique. At 
the end of the simulation, the full matter distribution in 
a light cone extending to z — 3 was saved in 417 time 
slices. Dark matter halos in this light cone were identi- 
fied with the Friends-of-Friends (FoF) algorithm using a 
linking parameter b — 0.2 (i.e. the linking length is one 
fifth of the mean interparticle separation of the simula- 
tion, or 195/i _1 kpc comoving). A 5 x 10 13 /i _1 Af Q cluster 
contains 744 particles. 

As shown in Figure [IJ the halo mass functions mea- 
sured from the simulated light con e agree well w i th the 
semi-analytic fitting formula from iJenkins et al.l (|2001l ) 
for a FoF linking length b = 0.2. The number of ha- 
los above a minimum mass per unit redshift per 100 
square degrees is plotted for three minimum masses 
M m j n = (1,2,4) x 10 14 M Q //i. The data points are mea- 
sured using all halos in the octant to minimize sample 
variance, but then are normalized to 100 square degrees. 
For any cluster survey, the minimum detectable mass is 
likely to be a function of redshift, and the three values 
chosen reflect the range expected for the upcoming SZ 
surveys. 



2.1.2. Adding Gas to Dark Matter Halos 

Gas is added to each clus ter following the p rescription 
of lOstriker et alj (|2005l) and lBode et~aTT(|2006l ). For each 
dark matter halo, the gravitational potential is found 
on a mesh with cell size 2e = 32.552/i~ 1 kpc. The gas 
density and pressure (or temperature) are found in each 
cell by assuming it is in hydrostatic equilibrium and has 
a polytropic equation of state with adiabatic index T — 
1.2. It is also assumed that initially the baryon fraction 
inside the virial radius is f2;,/f2 m , and that the energy per 
unit mass of the baryons equals that of the dark matter. 
To specify the central density and pressure with these 
assumptions requires two constraints. 

The first constraint is the pressure at the virial radius, 
which is the radius enclosing the virial overdensity cal- 
culated from spherical tophat co llapse (this is a slight 
change from lOstriker et al.l 120051 which used an over- 
density of 200 times critical). The dark matter kinetic 
energy is found in a buffer region extending for nine cells 
outside of the virial radius. This is then translated into a 
surface pressure, P s , by assuming the velocity dispersion 
is isotropic and the baryon fraction in this region equals 
the cosmic mean. The second constraint is conservation 
of energy. The final energy of the gas must equal this 
initial energy, after adjusting for any of the effects that 
can cause changes in energy. 

One such effect is that the gas distribution inside the 
cluster can expand or contract. This can, for example, 
reduce the baryon fraction inside the virial radius if the 
gas distribution expands further out. This would also 
reduce the energy of the gas as it does mechanical work 
pushing against the surface pressure. The surface pres- 
sure is assumed to remain at P s wherever the final gas 
radius ends up. 

Star formation will also affect gas energy. At z = Q, 
it is assumed that the stellar to gas mass ratio is 0.10. 
The star formation rate is a ssumed to follow a de layed 
exponential model (Eqn. 1 of lNagamine et aLll2006f ) with 
decay time t = 1 Gyr, so at higher redshift the star/gas 
mass ratio can be calculated. To make stars, the gas 
with the largest initial binding energy is removed, thus 
changing the total energy budget. Furthermore, some 
fraction of the rest mass turned into stars is taken to be 
converted into thermal energy in the gas, via supernovae 
and AGN. The energy added to the intracluster medium 
by these processes can be written as e/M*c 2 , where M* is 
the stellar mass and the feedback efficiency is estimated 
tobee/ = 5xl0~ 5 . This value is deter mined by fitting t o 
X-ray observations of nearby clusters (|Bode et alj r2006) ; 
at z=0, this added energy amounts to roughly 3 keV per 
particle. 

The upper panel of Figure O shows the resulting 
— T relation for all the clusters in the light cone 
at low redshift, z < 0.2. The temperature used is the 
X-ray emission- weighted temperature inside i?5oo- The 
kT > 5keV clusters are shown as open circles, and, for 
clarity, the median value of the clusters below 5 keV 
is shown as a line, with the shaded region enclosing 
90% of the clusters. Shown for comparison is data de- 
rived from X-ray clusters over th e same redshift range, 
taken from iReiprich fc Bohringerl ()2Q02f ) as adjusted by 
iMcCarthv et alj (|2004f ) . The simulated sample repro- 
duces the observed M 50 o — T relation quite well. The gas 
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Fig. 2. — The panels show the Afgno — T and Lx — T relations 
for the simulated clusters below z = 0.2 compared to those from 
observed X-ray clusters. The open circles are the simulated clusters 
above fcT ew = 5keV. The median value of the simulated clusters 
below 5 keV is represented by a solid line, with the shaded region 
enclosing 90% of the clusters. T ow is the emission- weighted X-ray 
temperature. 

temperature (or pressure) is fairly insensitive to the ex- 
act choice of e/ . Feedback has much more of an effect on 
gas density, making the Lx — T relation a more exacting 
comparison (|Bode et al.1 1200(1 ) . This relation is shown 
in the bottom panel of Figure O along with data po ints 
taken from t he ASCA cluster catalo g ((Horner 20011) as 
described in iMcCarthv et all (|2004D . A gain, the simu- 
lated clusters agree with the observed relation. 

2.2. Simulated SZ Signal 

For each of the simulated clusters, the gas prescription 
gives us the electron number density n e , temperature T e , 
and velocity v fields with which to model the SZ effect, 
including relativistic corrections. We assume that the 
gas is fully ionized with a helium mass fraction equal to 
0.24. The change in the microwave background intensity 
after passing through a path length dl in the direction h 
is given by 
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and the Y" 's and C"s are know n frequency-dependent co- 
efficients (|Nozawa et al.l (T998I) . The usual Compton y- 
parameter is given by 
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where dr = axn e dl is the optical depth through a path 
length dl. The first-order (non-relativistic) thermal and 
kinetic SZ signals are given by 
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respectively, where v\ os is the line-of-sight component of 
the peculiar velocity field. For a 10 keV cluster with 
v\ os — 1000 km/sec, the correction to the first-order 
thermal SZ from the 0(dr9 2 ) term at 145 GHz is about 
7.5%, but it becomes considerably more substantial near 
the null of the thermal SZ. The terms O(dTV\ QS e ) and 
O(dTV\ os 2 ) give abou t 8% and 1% corre ctions to the 
first-order kinetic SZ (|Nozawa et alj 119981 ). Note that 
the factor X 4 e x /(e x — l) 2 in equation [1] converts be- 
tween AT/Tcmb and AI u /Iq. 

Sky maps at the ACT observing frequencies of 145, 
215, and 280 GHz are made by tracing through the clus- 
ters and projecting them onto a cylindrical equal-area 
grid. Two strips of the sky, each 4 degrees wide in dec- 
lination and 360 degrees around in right ascension, are 
constructed and centered at 5 — — 5 and 8 — —55 degrees, 
respectively. Since the simulation light cone covers an oc- 
tant, only a quarter of the strip is unique and the other 
three-quarters are mirrored. The accompanying clus- 
ter catalog contains the following information: rcdshift, 
right ascension, declination, Mp F, ^200, -^500 , -R200, 
-R500, and integrated SZ and X-ray properties. (Mf f is 
the FoF derived mass and M200 and M500 are the cluster 
masses within -R200 and -R500 respectively, where the lat- 
ter are the radii at which the cluster mean density is 200 
and 500 times the critical density at the cluster redshift.) 
Figure [4] is a sample image showing the fully relativistic 
SZ signal at 145 GHz. This image is 33 degrees across 
in right ascension and 4 degrees wide in declination, cen- 
tered at 8 = —55 degrees. Black represents AT/Tcmb 
greater than —6 x 10~ 6 , and white represents AT/Tcmb 
less than —10 x 10~ 6 . The former value is chosen to be 
representative of ACT instrument noise in 0.2' pixels. 

In Figure [3l the first-order thermal SZ power spec- 
trum from the simulated map at 145 GHz is compared 
with a semi -analytic prediction derive d following the pre- 
scription of iKomatsu &; Seliakl ((2002). On large angular 
scales, / ~ 100, the two spectra are in agreement, but 
deviations on smaller scales are found. This is to be ex- 
pected, since our gas prescription accounts for star for- 
mation and feedback which reduce the gas density within 
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Fig. 3. — Solid curve is the first-order thermal SZ power spec- 
trum at 145 GHz from the simulation. Dashed curve is the first- 
order ther mal SZ power spe c trum at 145 GHz derived analyti- 
cally from Komatsu & Scliak (2002). For comparison, we include 
the lensed pri mary microwave backgrou nd power spectrum derived 
from CAMB l lChallinor &: Lewis 2005) for our input cosmology 
(dotted) and the power spectrum of infrared (dot-dashed) and radio 
(dash-dashed) point sources from the simulation, also at 145 GHz. 
The radio source power spectrum includes only sources fainter than 
35 mjy. Note that C; is in dimensionless units of (AT/T) 2 . 

clusters (jOstriker et al.ll200alBode et al.ll2006fl . resulting 
in a suppression of power on small (I > 1000) scales. In 
addition, the feedback pushes the gas out to larger radii 
and increases the overall temperature of the cluster, lead- 
ing to an enhancement of power on scales I < 1000. This 
demonstrates that one must have a good understanding 
of baryonic physics within clusters in order to extract cos- 
mology from the SZ power spectrum. Conversely, if the 
cosmology is well determined, then one can learn about 
cluster physics from the detected SZ signal. 

2.3. Simulated Primary CMB, Point Sources, and 
Galactic Dust 

A primary microwave background map is constructed 
to have the same large scale structure as we observe on 
the sky and small scale structure consistent with theo- 
retical expectation s. For our large scale map, we use the 
WMAP ILC map (jBennett et al.ll2003al iHinshaw et~aT1 
2006). For I < 20, the aim's are taken from the ILC map 
with no modification. At smaller scales, the ILC map 
is smoothed, so a Gaussian random field is added such 
that the ensemble average power equals the theoretical 
power spectrum taken from the WMAP LAMBDA web- 
site ( |http://lambda.gsfc.nasa.gov/product/map/dr2/| 
params / lcdm_all.cfmj^ We generate a map with 
HEALPix parameter N sidc = 4096 (totaling 201,326,592 
pixels). This is then interpolated onto the cylindrical 
coordinate system using the bilinear interpolation 
subroutine includ ed in the HEALPix distribution 
(jGorski et al.|[2005l ). 

The point sources are drawn from prescribed number 
counts. For a given flux bin of width AS, the number of 
sources in the survey is drawn from a Poisson distribution 



with mean \dN(> S)/dS\ASAQ. Within the bin, each 
source is assigned a random flux. We choose logarithmic 
bins that are narrow compared to the range of source 
fluxes. Each set of sources has flux counts tallied at some 
reference frequency vq. We scale the source flux from this 
frequency using a power law, S cx (f/fo) Q . Here a is a 
random variable, chosen for each source from a Gaussian 
distribution with mean a and variance tr^ un ique to each 
source population. For radio sources we use iKnox et al.l 
(|2004h number counts and choose a — —0.3 and a a = 0.3 
for scaling parameters. WMAP found a flatter spectrum 
with similar variance, but counts are expecte d to steepen 
above 100 GHz. For infrared sources we use lBorvs et al.l 
(2003) number counts and choose a — 3 and a a = 0.5. 

For the contribution of galactic ther mal dust emission, 
we us e the "model 8" prediction from iFinkbcinc r et all 
(1999), an ex trapolation to microw ave frequencies of the 
dust maps of lSchlegel et al.1 (|1998j ). This model is a two- 
component fit to IRAS, DIRBE, an d FIRAS data, and 
is shown by iBennett et al.l (|2003bf ) to be a reasonable 
template for dust emission in the WMAP maps. We 
create these galactic dust emission maps for completeness 
of the simulations. However, we do not currently include 
these dust maps when we investigate cluster detection. 

3. Y-M SCALING RELATION 

We would like to better understand the relation be- 
tween SZ flux and mass because this relation reflects clus- 
ter physics and also provides the link between SZ clus- 
ter catalogs and cosmological parameters, which are con- 
strained by the cluster distribution as a function of mass. 
A more informed understanding about how changes in 
cluster physics and cosmology alter the Y — M relation 
will allow for more information about cosmological pa- 
rameters to be derived from cluster surveys. 

To investigate the SZ flux - mass cluster scaling rela- 
tion, the integrated Compton-y parameter is compared 
to cluster mass for the clusters in our simulation. The 
integrated Compton-y parameter is the Compton-y pa- 
rameter integrated over the face of the cluster, i.e. Y = 
(J ydfl)d A (z), where d 2 A (z) is the angular diameter dis- 
tance. In the self-similar model, a virialized halo of mass 
M v i r has a virial temperature equal to 

T vir cx [M vir £(z)] 2 / 3 , (8) 

where for a flat ACDM cosmology 

E{z) = [tt m {l + zf + n K ] l X. (9) 

If clusters were isothermal, we would expect them to sat- 
isfy the relation 

Y cx / gas M halo T, (10) 
where / gas is the cluster gas mass fraction. Thus we find 

YKf s&s M r tE{zf'\ (11) 

for the self-similar SZ flux - mass scaling relation. Since 
clusters are not isothermal and not always in virial equi- 
librium, deviations from the self-similar relation are ex- 
pected. Moreover, variation of / gas with redshift and 
cluster processes such as star formation and feedback can 
also cause deviations from self-similarity. 

From the cluster catalog described in §2, we choose 
a sub-sample with M200 > 7.5 x 10 13 Af Q to stay well 
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Fig. 4. — A sample image showing the fully relativistic SZ signal at 145 GHz. The image dimensions are 33 degrees across in right 
ascension by 4 degrees wide in declination centered at S = —55 degrees. Black represents AT/Tcmb greater than —6 X 10~ 6 , and white 
represents AT /Tqmb l ess than —10 X 10~ 6 . The former value is chosen to be representative of ACT instrument noise in 0.2' pixels. 



above our catalog mass limit of 5 x 10 Mq and ensure 
catalog completeness. This leaves us with a sample of 
about 10 5 clusters. For each cluster, we calculate Y200, 
which is the projected SZ Compton-y parameter in a 
disk of radius i?200- Clusters are then ranked by mass, 
and for bins of 250 clusters each, the mean of M200 as 
well as the mean and standard deviation on the mean of 
Y200 x E(z)~ 2 / 3 are calculated in each bin. These values 
are plotted in Figure [5] We then fit these values to the 
power-law relation 

^200 _ 1 n/3 ( -^200 
Vl0 14 M ra 



(12) 



We find a = 1.876±0.005 and /? = -5.4774±0.0009 with 
a reduced \ 2 of 1.004. This suggests that the Y — M re- 
lation is close to a power-law. This slope is steeper than 
that for the self-similar model given in equation [TT] This 
steepening is expected since we assume feedback is inde- 
pendent of total cluster mass. Energy input of roughly 3 
keV per particle will clearly have a relatively larger effect 
on small clusters with lower virial velocities than on more 
massive ones. The feedback reduces the gas mass frac- 
tion more in lower mass clusters, as more gas is pushed 
out in proportion to the cluster's total mass (i.e., the gas 
fraction inside i?2oo is smaller). As a result, Y200 is also 
decreased more in the lower mass clusters. The power- 
law index that we find is slightly steeper t han that foun d 
by some hydrody n amic simulat i ons (e .g. iNag ai (2006); 
iMotl et all (120051): IWhite et all (|2002D ). This is under- 
standable since we include feedback from active galac- 
tic nuclei as well as supernovae, which is different from 
the hydrodynamic simulations, and the increased level of 
feedback steepens the slope. It has been estimated that 
the feedback from active galactic nuclei may be roughly 
an order of magnitud e larger than that from supernovae 
(|Ostriker et al.ll2005l ). so it potentially has a significant 
eff ect. The slope we fin d is also in between that found 
by Ida Silva et al. (2004) for their hydrodynamic simula- 
tions including radiative cooling alone and cooling plus 
preheating. 

A departure from a power-law is observed when we re- 
strict our cluster sample to have higher minimum masses. 
As the minimum mass of the sample is increased, the 
power-law index gets flatter, illustrating the important 
effect of feedback on low-mass clusters as well as the fact 
that the lowest mass clusters dominate the slope of the 
Y — M relation because their numbers increase rapidly as 
mass is decreased. We illustrate this departure in Table 
2 where we have preformed a similar analysis as above 
for clusters with M200 > 2 x 10 14 M Q , this time using 
bins of 20 clusters each. We find a = 1.81 ± 0.02 and 
= -5.463 ± 0.009 with a reduced X 2 of 1.14. This de- 
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Fig. 5. — Y200 X E(z)~' 2 l' i vs. A/200 for » 10 5 clusters within the 
N-body plus gas simulation. Each point represents the mean value 
for 250 clusters, and error bars represent the error on the mean. 



The parameters and reduced 
Table [I] 



X for the best-fit line are given in 



parture indicates that while the Y — M relation is close 
to a power-law, there is some curvature. 

We also divide our clusters above M200 > 0.75 x 
10 14 M Q into four different redshift bins (z < 0.3, z G 
(0.3, 0.6), z G (0.6,0.9), and z G (0.9,1.5)) to see the 
evolution of the power-law index and normalization with 
redshift. The best-fit lines to Y200 x E(z)~ 2 / 3 versus 
M200 are shown in Figure [H and the best fit a and (3 are 
given in Table [T] Some slight evolution of the power-law 
index and normalization with redshift is apparent. This 
is expected since clusters are not perfectly self-similar. 
Furthermore, the stellar to gas mass ratio is a function of 
redshift in our simulations, and thus the feedback, which 
is proportional to the star formation, is also a function 
of redshift. Since more feedback lowers Y disproportion- 
ately at the low-mass end, we find the steepest slope for 
the clusters with z < 0.3, which have undergone more 
star formation and feedback. We find the flattest slope 
for the clusters in the highest redshift bin. Examining 
only the higher mass clusters, M 20 o > 2 x 1O 14 M , we 
see little obvious evolution with redshift of the slope and 
normalization considering their respective error bars (cf. 
Table 2) . This again suggests that the higher mass clus- 
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Fig. 6. — Y200 X E(z) -2 / 3 vs. M200 for four different redshift 
regions. As in Figure 5, each point represents the mean of 250 
clusters, and error bars represent the error on the mean. The 
parameters and reduced x 2 f° r the best-fit lines are given in Table 

m 

ters are less sensitive to feedback processes. 

The scatter that we find for this Y — M relation is 
most likely underestimated, as we assume that the gas is 
described by hydrostatic equilibrium, a single constant 
polytropic index, and the same amount of feedback per 
stellar mass for all clusters. Also, we have not included 
other sources of non-thermal pressure support which may 
contribute non-negligibly to the total press ure. Some of 
these effects, such as turbulent pressure (see lRasia et al.l 
(2004)), are included in all current high-resolution hy- 
drodynamic simulations of galaxy clusters. Initial at- 
tempts have also b een made to includ e such effects as 
magnetic fields (e.g. iDolag et al.l (|1999t )1 and cosmic rays 
(|Pfrommer et al.ll2006f K Hvdrodvnamic simulations can 
provide a robust check and calibration for semi-analytic 
models, but continuing work is required to include all 
of the relevant gas physics in both methods. Possible 
changes in the parameters of the predicted Y — M rela- 
tion can be expected as models achieve greater realism. 

4. CLUSTER DETECTION 

To study cluster detection, three maps of the ACT 
strip are first created (combining the primary microwave 
background, infrared and radio point sources, and the 
full SZ effect) at the frequency channels of 145, 215, and 
280 GHz. The ACT instrument is simulated by convolv- 
ing the ACT strip maps at 145, 215, and 280 GHz with 
a Gaussian beam of full- width at half-maximum equal to 
1.7', 1.1', and 0.93' respectively. Gaussian random noise 
is then added to the maps of a equal to 2/iK, 3.3/iK, and 
4.7/iK per beam for the three frequencies respectively. 
The above are preliminary detector noise estimates for 
extensive dedicated observations of the ACT strip. We 
produce a single map of the SZ clusters by applying a 
multi- frequency Wiener filter to the three ACT maps. 



This filter requires a model power spectrum of the first- 
order thermal SZ signal from the clusters (see Eqn. [6J. 
We first use a spectrum derived from the simulations 
themselves, and t hen one derived from the a nalytic pre- 
scription given bv lKomatsu fc Seljakl (|2002l ) to test the 
sensitivity of the filter performance on the input clus- 
ter model. Clusters are then identified within the fil- 
tered map using a simple peak-finding algorithm. These 
clusters are matched to a catalog documenting the in- 
put map. From this matching process, we determine the 
completeness and purity of our detected cluster sample. 

4.1. Multi- Frequency Wiener Filter 

In the literature, a number of different filters have been 
employed to extract the SZ signal from microwave simu- 
lations (see introduction for references). Two commonly 
employed filters are matched or scale-adaptive filters and 
Wiener filters. Matched filters make an explicit assump- 
tion about cluster profiles while Wiener filters assume 
knowledge of the SZ power spectrum. Both have qual- 
itatively similar shapes: they suppress the primary mi- 
crowave background at large scales and noise at small 
scales. Our motivation for choosing to use a Wiener fil- 
ter stems from the fact that it recovers the minimum 
variance SZ map and is simple to imple ment. We com- 
pare our results with iMelin et alj (|2006l ) who employ a 
scale-adaptive filter and find a Wiener filter yields com- 
parable results. We also use two different Wiener filters 
and find the differences are minor. 

Below we describe the map filter in detail. We wish to 
recover one filtered map (of clusters) from several ACT 
maps (at three different frequencies). The latter are de- 
scribed by the data vector d. We write d = Rs + n. 
where R is the instrument response to the Compton-y 
signal map s (see Eqn. [5]), and n is the noise. In this 
notation, R includes both the frequency dependence of 
the first-order thermal SZ signal (see Eqn. [6]) and the 
convolution with the beam. The noise n includes detec- 
tor noise, primary microwave background, point sources, 
and contributions to the SZ signal other than from the 
first-order thermal SZ. We denote the covariance of the 
signal, noise, and data by S = (ss^), N = (nn^), and 
D = (cUf) = RSR^ + N, assuming signal and noise are 
largely uncorrelated. 

The multi-frequency Wiener filter is given by 



W = SR f [RSR f + N] 1 = SR f D- 



(13) 



and it gives the least-squares signal map reconstruction 
via s = Wd. Here we are assuming statistically homoge- 
neous signal and noise, so all the covariances are (block) 
diagonal in the Fourier representation: they are the auto 
and (frequency) cross-spectra. This makes the appli- 
cation of the Wiener filter straightforward: we Fourier 
transform the ACT maps, apply the filter, and inverse 
Fourier transform to obtain the SZ cluster map. (Sec 
iTegmark fc Efstathioul (|1996[ ) for a useful reference on 
Wiener filtering.) 

4.2. Filtering Maps 

To filter the maps, we divide the ACT strip (with SZ, 
primary microwave background, and point sources com- 
bined) into roughly 4° x 3.25° patches. The above dimen- 
sions are chosen to give square patches in pixel space. 
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TABLE 1 

Best-fit parameters for the best-fit lines depicted in Figures [5] and [6] fit using the power-law given in equation [12] 
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TABLE 2 

Best-fit parameters for clusters with M200 > 2 x 10 14 M@, fit using the power-law given in equation [T2l 



To compute D, the power and cross-power spectra for 
each patch is calculated and averaged over the whole 
strip. For the first model cluster power spectrum (S), we 
compute the spectrum from the simulations in a similar 
fashion. Later we also us e the semi-analytic spectrum of 
IKomatsu fe Sella! (|200l . 

To avoid aliasing of signal at the edges of each patch, 
we use an overlap-and-save method to filter overlapping 
pieces of the map. This permits us to discard all the 
vertical edges (edges of constan t right ascension) on the 
patch borders (|Press et al.| [l997). In this way, we remove 
any discontinuities in the filtered map at the boundaries 
of the patches. The horizontal edges are less of a con- 
cern since they can be discarded after the whole strip is 
filtered. This filtering procedure provides a filtered map 
of the entire ACT strip consisting of the recovered SZ 
signal. 

4.3. Identifying Clusters 

Clusters are identified by searching for all peaks in the 
filtered map above a given threshold value. The thresh- 
old cuts we apply are 1, 2, and 3 times the standard 
deviation of the filtered map. A peak is identified simply 
as a pixel larger than any pixel within a radius r = 6p 
around it, where p is the pixel size (p ~ 0.21' for these 
maps). We choose this radius both because w 1.3' is a 
typical cluster size and in an effort to be consistent with 
iMelin et all (pOOl for comparison purposes. Once a list 
of peaks is compiled, the peaks are matched to clusters 
in the input cluster catalog. We identify a match if the 
peak and catalog cluster are within a distance r of one 
another. Multiple peaks matching a single catalog clus- 
ter are allowed, as are multiple catalog clusters matching 
a single peak. Follow-up observations to determine clus- 
ter redshifts should sort these cases out. Any peak that 
is not within r of any catalog cluster, is flagged as a false 
detection. 

It is possible that some peaks match catalog clusters 
just by chance alignment. However, since clusters down 
to only 5 x 10 13 M Q are included, all the clusters in the 
simulation are comparable to or above ACT instrument 
noise, and thus are in principle detectable. So we accept 



that some small fraction of the simulation clusters may 
be detected by chance alignment with spurious peaks and 
do not try to distinguish these cases from real matches. 

4.4. Results 

First we examine the case when the three-frequency 
ACT strip maps contain no radio or infrared point 
sources. The completeness and number of detected clus- 
ters per deg 2 arc given in Figure [7] as a function of both 
Af2oo and I200 (the latter measured as a solid angle in 
units of arcmin 2 ). The dashed, dotted, and dash-dotted 
curves represent the 1-cr, 2-er, and 3-cr threshold cuts re- 
spectively, and the solid line represents the simulation 
mass function. The purity of each cluster sample is given 
in the legend (purity = 1 — % of false-positive detections; 
completeness = number of detected clusters/number of 
input clusters). In the absence of point sources, using 
a 3-cr threshold cut, our sample of detected clusters is 
96% complete down to 2 x 10 14 M Q with only 4% con- 
tamination from false positives. The completeness as a 
fu nction of Y200 and p urity are comparable to that found 
bv lMelin et al.1 (|2006[ ) when they investigated a no-point- 
source case for the South Pole Telescope (SPT) using 
a matched filter and simulated clusters that perfectly 
matched their filter. 

In the completeness versus I200 plot, there is a sharp 
decrease in the completeness from 1 to 0.8 at Y200 = 
2 x 10 3 arcmin 2 . This effect is caused by a single cluster 
at z « 0.04 with A/ 20 o ~ 8 x 1O 13 M . Since this cluster 
is so nearby, it extends over many pixels, and thus has a 
large Y200 i n units of arcmin 2 . However, its signal in each 
pixel is below that of the 1-ct threshold cut, so it is not 
detected. If the plot had been made with Y200 in units 
of physical area (Mpc 2 ), this feature would disappear. 
Therefore, this feature is an artifact reflecting the fact 
that arcmin 2 is not a unit intrinsic to a cluster. It also 
suggests that for detecting extended, low-redshift clus- 
ters, searching for a peak pixel is not the best method. 
Profile-fitting may provide a better alternative in this 
case. 

Next, infrared point sources are added to the three- 
frequency ACT maps following the prescription de- 
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Fig. 7. — Completeness and purity of the detected cluster sample with ACT instrument specifications, in the case of no point sources. 
Different curves represent 1-tr, 2-<r, and 3-cr threshold cuts in the detection algorithm, where a is the standard deviation of the filtered 
map. The solid curve represents the simulation mass function. The purity of each cluster sample is given in the legend. See text for further 
details. 



scribed in §2.3. Infrared sources up to a flux limit of ~ 0.4 
Jy at 145 GHz are included since the number of sources 
with higher flux is effectively zero for the area of the 
ACT strip given the above prescription. Figure [8] shows 
the cluster detection results after adding this population 
of point sources. For this case, using a 3-cr threshold cut, 
the detected cluster sample is 95% complete down to 
3 x 10 14 M Q with 14% contamination from false-positive 
detections. The perc entage of fal s e-posi tives is compara- 
ble to that found bv lMelin et alj (|2006l ) when they used 
N-body, as opposed to spherical, clusters for their no- 
point-source case. A direct comparison of completeness, 
including infrared point sources, with their work is diffi- 
cult because we do not use isothermal, spherical clusters 
and we include clusters down to a lower mass limit, which 
crowds the field and increases confusion. 

We further investigate the case where both infrared 
and radio point sources are added to the three-frequency 
ACT strip maps. The radio sources follow the prescrip- 



tion described in §2.3. Figure [5] shows the cluster de- 
tection results including radio sources in the ACT maps 
below three different flux limits (3.5 mJy, 35 mJy, and 
350 mJy at 145 GHz). All the curves correspond to 
cluster detection using a 3-cr threshold cut. The dashed 
curve shows that if all the radio sources above 3.5 mJy 
at 145 GHz are removed from the three-frequency ACT 
maps, the infrared-point-source-only case is recovered, 
depicted in Figure [H If all the radio sources above 35 
mJy at 145 GHz are removed from the three ACT maps, 
we obtain results that are of only slightly lower quality 
than the previous case. These results can be understood 
by noting that the inclusion of radio sources below 35 
mJy and infrared sources increases the average power in 
the simulated ACT strip on scales of a few arcminutes 
by roughly an order of magnitude, as compared to not 
including point sources. The inclusion of all the radio 
sources above 35 mJy increases the power on these scales 
by several orders of magnitude more. This large increase 
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Completeness and purity of the detected cluster sample in the case of infrared point sources only. 



in noise power is why these few very bright sources limit 
cluster detection with this technique. For an arcmin 2 
beam at 145 GHz, 3.5 mJy and 35 mJy correspond to 
w 100/^K and « 1000/zK respectively. Thus achieving the 
removal of radio sources above the former flux limit may 
require interferometric observations of the ACT strip at a 
lower frequency (e.g. with the Australia Telescope Com- 
pact Array (ATCA) at w 20 GHz). Achieving the latter 
flux limit of radio sources is most likely possible with 
ACT alone. (Sources of w 1000//K should be easy to 
identify, and we find 700 of them in our simulations 
of the ACT strip, which is not an unmanageable number 
to remove.) Assuming removal of radio sources above 
35 mJy at 145 GHz, our detected cluster sample is 90% 
complete down to 3 x 10 14 M Q and 84% pure. 

The filter used in obtaining the above results was cre- 
ated using the Compton-y power spectrum from the sim- 
ulations themselves. Since we shall not know beforehand 
the true cluster power spectrum, we remake the filter 
using a Compton-y power spectrum derived analytically 



from the prescription in lKomatsu fc Seliakl ((2002). With 
this filter, we repeat the above exercise. Figure [10] dis- 
plays the cluster detection results using this filter and 
adding both infrared and radio point sources. The dif- 
ferent curves again correspond to the three different flux 
cuts of radio sources, and a 3-cr threshold cut is used for 
them all. For the three flux cut cases, there is a slight in- 
crease in completeness of the cluster sample as compared 
to Figure [9] However, there is also a slight decrease in 
the purity of the cluster sample. If we compare the first- 
order thermal SZ power spectrum from the simulations 
and derived analytically (Fig. [3]), we see that the fil- 
ter made with the analytic power spectrum leaves more 
small scale power in the filtered map. This increases the 
sample completeness, but also decreases its purity, as 
more point sources are left in the map as well. Overall, 
however, the difference in using the analytic power spec- 
trum in the filter, as compared to the power spectrum 
from the simulations, is small. This suggests that the 
performance of the filter is not overly sensitive to how 
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the cluster physics is modeled when constructing it. 
5. FUTURE WORK 

In this work, we did not correlate the radio and in- 
frared point sources with the clusters. This was due 
to a lack of observational data regarding these correla- 
tions. If either population is both strongly correlated 
with clusters and exhibits a relatively flat spectrum be- 
tween the infrared/radio regime and the microwave, then 
some further percentage of clusters may fail to be de- 
tected. Recent BIMA and OVRO observations targeting 
massive galaxy clusters at 28.5 GHz suggest an overabun- 
dance of mJy radio sources in these massive clu sters as 
opposed to non-cluster fields (jCoble et al.ll2006l ). How- 
ever, they also find a steeper spectral index for these 
sources as compared to the brigh ter, rarer popu l ation o f 
radio sources found by WMAP (jBennett et alj l2003bf ) . 
If their spectral index steepens further above 30 GHz, 
these former sources may not be an issue. But if their 



slope flattens out, some of these sources may fill in the 
SZ decrement of some percentage of clusters at 145 GHz. 
Without a better understanding of these correlations and 
observations of these point sources in the range of 145 to 
280 GHz, it is difficult to model both infrared and radio 
sources more realistically. Upcoming microwave surveys 
and follow-up or concurrent radio and infrared observa- 
tions overlapping the survey areas will hopefully clarify 
these current unknowns. It is useful to keep in mind 
that, for constraining cosmology, obtaining a large num- 
ber of detected clusters is less important than obtaining a 
solid understanding of the selection criteria. So if point 
sources affect cluster detection more significantly than 
anticipated, as long as this effect can be accurately char- 
acterized, it should not pose an overwhelming hurdle for 
cosmological investigations. 

Some further avenues of investigation that we did not 
cover in this work include the following. The instrument 
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Fig. 10. — Completeness and purity of the detected cluster sample for the cases described in Figure [9] but this time using a filter 
constructed with a cluster Compton-y power spectrum derived analytically from Komatsu & Scliak (2002), instead of taken directly from 
the simulations. The completeness is slightly increased and the purity slightly decreased as compared to the previous Figure, but overall 
the difference is small. 



noise can be more realistically modeled than we have 
done here, and atmospheric noise can be included. The 
lensing of the primary microwave background and point 
sources by clusters can be modeled to study its effect on 
cluster detection and SZ flux recovery. Moreover, the 
optimal recovery of cluster SZ flux, which is an issue dis- 
tinct from cluster detection, should be investigated. This 
will aid in flux - mass identifications and SZ component 
separation. Some r ecent work addres sing the later is- 
sue can be found in iMelin et~aT1 (|2006t ) . We also by no 
means suggest that Wiener filtering is the ideal cluster 
detection tool, and further study in this direction is also 
warranted. 

6. CONCLUSION 

New microwave instruments have the potential to de- 
tect many galaxy clusters out to redshift z=l and be- 
yond by their SZ signature. These galaxy clusters pro- 
vide information about the growth of structure, which in 



turn constrains cosmology. However, the SZ signal is em- 
bedded within the signals from the primary microwave 
background, radio and infrared point sources that still 
have significant flux in the microwave regime, and galac- 
tic emission. Thus, in an effort to study the detection of 
clusters via their SZ effect, we simulated the microwave- 
sky over both the proposed ACT observing region and 
a parallel strip centered at S = —5 degrees. To realisti- 
cally model the small scale cluster physics, which affects 
the SZ signal, we comb ined an N-body simulation with 
an an alytic gas model (|Ostriker et al.ll2005t iBode et all 
2006). This allowed us to model cluster processes such 
as star formation and feedback over much larger volumes 
than can currently be achieved by hydrodynamic simu- 
lations. One product of these cluster simulations is that 
we are also able to study features in the SZ flux - mass 
scaling relation, such as its slope and normalization and 
their evolution with redshift. Knowledge of this relation 
is necessary to connect SZ flux to cluster mass, which is 
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the quantity of interest for cosmology. 

In investigating the relation between SZ flux and mass, 
roughly 10 5 clusters from our N-body plus gas simula- 
tions are fit to the power-law relation given in equa- 
tion [l2l The best-fit power-law index is found to be 
a = 1.876 ± 0.005, and the reduced \ 2 to be 1.004. This 
slope is steeper than 5/3, which is the expectation for 
the self-similar model. This steeper slope is consistent 
with the inclusion of feedback from both active galactic 
nuclei and supernovae in these simulations which serves 
to lower the SZ flux more for lower-mass clusters. This 
slope is understandably steeper than that for hydrody- 
namic simulations, which do not include feedback from 
active galactic nuclei and thus have a lower amount of 
overall feedback than in this simulation. We find some 
redshift dependence for both the power-law index and 
normalization of this relation. This is expected since 
clusters are not in fact self-similar, and because star for- 
mation and cluster feedback are redshift dependent in 
our simulations. The steepening of the power-law index 
at lower redshifts is consistent with more star formation 
and feedback occurring at later times. Considering only 
the higher mass clusters, M200 > 2 x 10 14 M©, we find 
a flatter slope of a = 1.81 ± 0.02 with a reduced % 2 of 
1.14. We also see less obvious evolution of the slope and 
normalization with redshift. This is consistent with the 
higher mass clusters being less sensitive to feedback pro- 
cesses. It should be noted that the details of this Y — M 
relation can change as additional sources of non-thermal 
pressure support such as magnetic fields, cosmic rays, 
and turbulence are included in semi-analytic gas models. 

Our projections for cluster detection with the ACT 
instrument arc also promising. Cluster detection was in- 
vestigated under varying levels of point source contami- 
nation utilizing a multi-frequency Wiener filter and peak- 
finding algorithm. These results suggest that in the ab- 
sence of point sources, considering only peaks above 3-cr 
in the filtered map, ACT can obtain a cluster catalog that 
is 96% complete above 2 x 10 14 M Q and 96% pure (4% 
false-detections). The inclusion of infrared sources re- 
sults in a catalog that is 95% complete above 3 x 10 14 M Q , 
and 86% pure. When all radio sources are included in ad- 
dition, with no flux cutoff, we find that the noise from the 
brightest radio sources interferes significantly with clus- 
ter detection using this technique. However, if the bright- 



est sources are removed from the ACT maps, there is a 
considerable improvement. Removing all radio sources 
above 35 mJy at 145 GHz from the three-frequency ACT 
maps results in a catalog that is 90% complete above 
3 x 10 14 M© and 84% pure. Removing all radio sources 
above 3.5 mJy at 145 GHz gives the same results as the 
inclusion of only infrared point sources. These results 
are encouraging, since 35 mJy sources should be easily 
identifiable and removable from ACT maps without need 
for further observations at alternate frequencies. 

This study was repeated using a cluster model in our 
filter that differed from that in our simulations. Doing 
this produced only a small change in the cluster detection 
results. This suggests that the filter performance is not 
very sensitive to the details of the cluster model used to 
create it. As a result, efforts to characterize the selection 
function are made easier. 

The potential of new microwave instruments to pro- 
vide knowledge about cluster physics and cosmology is 
substantial, as our results further confirm. It is hoped 
that the simulations created to pursue the above studies 
will serve as a useful tool for future investigations toward 
making this potential a reality. 
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